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Recent neutron scattering experiments on aqueous salt solutions of amphiphilic t-butanol by 
Bowron and Finney [Phys. Rev. Lett. 89, 215508 (2002); J. Chem. Phys. 118, 8357 (2003)] 
suggest the formation of t-butanol pairs, bridged by a chloride ion via O— H ■ • ■ Cl~ hydrogen-bonds, 
and leading to a reduced number of intermolecular hydrophobic butanol-butanol contacts. Here we 
present a joint experimental/theoretical study on the same system, using a combination of molec- 
ular dynamics simulations and nuclear magnetic relaxation measurements. Both MD simulation 
and experiment clearly support the more classical scenario of an enhanced number of hydrophobic 
contacts in the presence of the salt, as it would be expected for purely hydrophobic solutes [J. Phys. 
Chem. B 107, 612 (2003)]. Although our conclusions arrive at a structurally completely distinct 
scenaxio, the molecular dynamics simulation results are within the experimental errorbars of the 
Bowron and Finney work. 



I. INTRODUCTION 

Nonpolar solutes, such as noble gases or alkanes, don't 
like to be dissolved in water. Consequently, they are con- 
sidered as "hydrophobic" and their corresponding solva- 
tion free energy is found to be large and positive [1-5]. 
This effect is typically found to be significantly strength- 

I ened when salt is added, leading to a further reduced 
solubility of hydrophobic species such as noble gases, or 
methane [6, 7]. The increasing excess chemical potential 
is usually found to be proportional to the salt concen- 
tration over large concentration ranges and is therefore 
parameterized in terms of Setschenow's concentration in- 
dependent salting out coefficient [6]. In line with the 

' observation of an increased positive solvation free energy 
upon addition of salt, Ghosh et al. [8] report an increased 
number of hydrophobic contacts in a diluted aqueous so- 
lution of methane. Moreover, for the case of hydrophobic 
interactions in a hydrophobic polymer chain, an approxi- 
mately linear relationship between the salt concentration 
and the strength of pair-wise hydrophobic interactions 

] has been determined [9]. This observation seems to be 
in line with the general finding of Koga et al., that the 

] excess chemical potential of a small hydrophobic particle 
and the strength of hydrophobic pair interactions appears 
to be (almost) linearly related [5, 10]. 

However, purely hydrophobic compounds are proba- 
bly not very typical representatives of biophysical con- 
stituents. Usually, proteins and membrane-forming lipids 
are amphiphilic in the sense that they are composed of 
both, hydrophobic and hydrophilic groups, where the lat- 
ter ones ensure a sufficiently high solubility in an aque- 
ous environment. In addition, the delicate interplay be- 
tween hydrophobicity and hydrophilicity is exploited by 
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nature to control the dimensions of molecular aggregates 
in aqueous solution [11-14]. Small amphiphilic alcohols 
might thus be considered as a minimalist model system 
to explore the subtle interplay between hydrophobic and 
hydrophilic effects. Numerous experimental studies on 
alcohol aggregation in aqueous solution have been re- 
ported, based on nuclear magnetic resonance [15-19], as 
well as light-, x-ray-, and neutron scattering techniques 
[20-26]. 

Salts are known to influence a number of properties of 
aqueous solutions in a systematic way [27] . The effect of 
different anions and cations appears to be ordered in a 
sequence, already proposed by Hofmeister in 1888 [28], 
deduced from a series of experiments on the salts abil- 
ity to precipitate "hen-egg white protein" . However, the 
exact reason for the observed specific cation and anion 
sequences is still not completely understood, since the 
same salt that can precipitate a protein at one concen- 
tration can 'salt it in' at another [29]. Model calculations 
[30], as well as nuclear magnetic relaxation experiments 
[16] propose a delicate balance between ion adsorption 
and exclusion at the solute interface, tuned by the sol- 
vent (water) structure modification according to the ion 
hydration [31, 32] and hence possibly subject to molecu- 
lar details. 

Recently, Bowron and Finney provided a detailed 
mechanistic picture of the possible salting out process 
of t-butanol (TEA) in aqueous solution. The atomistic 
structure of the solution was determined from neutron 
scattering experiments [25, 26], varying solute and sol- 
vent isotopic compositions [33]. However, their analysis 
relies largely on the accuracy of the employed empiri- 
cal potential structure refinement (EPSR) technique of 
Soper [34, 35]. Their main observation is that TBA- 
molecules form dimers that are connected by hydrogen 
bonds to a central chloride anion. Salting out appears 
hence due to solute aggregates, which are formed by 
anion-bridges between the hydroxyl-groups, increasing 
the solutes overall hydrophobic surface and thus reducing 
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the solubility of the whole complex [36-38]. A straight- 
forward conjecture would suggest that the salting out of 
proteins could be driven by analogous anion-bridged ag- 
gregates. We would like to point out that the proposed 
mechanism has remarkable similarity with the "differen- 
tial hydrophobicity" concept of Burke et al. [39] used 
to explain the specific protein-aggregation behavior ob- 
served in Huntington's disease. 

The molecular dynamics simulations discussed here, 
however, do not show any evidence for a "salting out"- 
scenario as proposed by Bowron and Finney. Instead, 
upon addition of salt we find an increased number of hy- 
drophobic contacts of the TBA molecules which increases 
with higher salt concentration. Using a combination of 
molecular dynamics simulations and nuclear magnetic re- 
laxation experiments we show that an association param- 
eter based on NMR measurable quantities and introduced 
by H.G. Hertz et al. [40, 41] is a useful measure for the 
association of TBA molecules in the present case. Both 
simulation and NMR experiments consistently support 
the classical picture of an enhancing hydrophobic associ- 
ation in the case of aqueous TBA/salt solutions. 



II. METHODS 

A. Dipolar nuclear magnetic relaxation and 
correlations in the structure and dynamics of 
aqueous solutions 

The molecular dynamics simulations yield the time- 
dependent positions of the atomic nuclei. Experimen- 
tally, the individual molecular positions are not avail- 
able. However, a useful measure of molecular association 
is experimentally accessible via the measurement of nu- 
clear spin relaxation rates and is discussed in a separate 
section below. In this section we would like to briefly 
summarize the underlying theory. 

The most important contribution to the relaxation rate 
of nuclear spins with / = 1/2 is the magnetic dipolc- 
dipole interaction. The relaxation rate, i.e. the rate at 
which the nuclear spin system approaches thermal equi- 
librium, is determined by the time dependence of the 
magnetic dipole-dipole coupling. For like spins, it is [42] 

Tfi = 274^^/(7 + l)(/xo/47r)2 (1) 

^ A f \^ Do,2[^ij{0)] £>0,2[%(t)]\ 

where D}.^m[^] is the fc,m-Wigner rotation matrix ele- 
ment of rank 2. The Eulerian angles 0(0) and Vl{t) at 

time zero and time t specify the dipole-dipole vector rel- 
ative to the laboratory fixed frame of a pair of spins, 



denotes their separation distance and /ig is the permit- 
tivity of free space. The sum indicates summation of all 
j interacting like spins in the entire system. (. . .) denotes 
averaging over all equivalent nuclei i and all time zeros. 
For the case of an isotropic fluid and in the extreme nar- 
rowing limit Eq. (1) simplifies to [43] 

00 

Tf 1 = 2^^h^I{I + 1) (^) Y G2(i) dt . (2) 



The dipole-dipole correlation function here is abbreviated 
as G2{t) and is available through [43, 44] 

G2{t) = ^E^^'(0)^fW^2[cos%(t)]^ , (3) 

where dij (t) is the angle between the vectors joining 
spins i and j at time and at time t [43] and P2 is the 
second Legendre polynomial. 

To calculate the integral, the correlation function G2 {t) 
can be separated into an r~^-prefactor, which is sensi- 
tive to the structure of the liquid (average internuclear 
distances) and a correlation time T2, which is obtained 
as the time-integral of the normalized correlation func- 
tion G2{t), and which is sensitive to the mobility of the 
molecules in the liquid, 

jG2it)dt = . (4) 

From MD-simulation trajectory data the correlation 
function G2 and hcmcc^ Ti can be calculated directly. 
From the definition of the dipole-dipole correlation func- 
tion in Eq. (3) it follows that the relaxation time Ti is 
affected by both, reorientational and translational mo- 
tions in the liquid. Moreover, it is obvious that it also 
depends strongly on the average distance between the 
spins and is hence sensitive to changing inter- and in- 
tramolecular pair distribution functions [40, 41]. In addi- 
tion, the r~^-wcighting introduces a particular sensitivity 
to changes occurring at short distances. For convenience, 
one may divide the spins j into different classes accord- 
ing to whether they belong to the same molecule as spin 
i, or not, thus arriving at an inter- and intramolecular 
contribution to the relaxation rate 

^1 ^ = ^1, inter + ^1, intra' C^) 

which are determined by corresponding intra- and in- 
termolecular dipole-dipole correlation functions 6*2, intra 
and G2,inter- The intramolecular contribution is basi- 
cally due to molecular reorientations and conformational 
changes and has been used extensively to study the re- 
orientational motions, such as that of the H-H- vector in 
CHa-groups in molecular liquids and crystals [45]. In 
the course of this paper, however, we are particularly in- 
terested in the association of solute molecules, and will 
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therefore focus on the intermolecular contribution (see 
also experimental section). 

The structure of the liquid can be expressed in terms 
of the intermolecular site-site pair correlation function 
gij{r), describing the probability of finding an atom of 
type j in a distance r from a reference site of type i 
according to [46] 



1 



NiPj 



(6) 



\fe=i ;=i 



where pj is the number density of atoms of type j. The 

prcfactor of the intermolecular dipolc-dipolc correlation 
function is hence related to the pair distribution function 
via an integral of the pair correlation function 



E^^'(O)) = j 9ij{r) 47rr2 dr 



(7) 



Since the process of enhanced association in a molecu- 
lar solution is equivalent with an increase of the nearest 
neighbor peak in the radial distribution function, Eq. 7 
establishes a quantitative relationship between the de- 
gree of intermolecular association and the intermolecular 
dipolar nuclear magnetic relaxation rate. 



B. Self-association: the Hertz "A" -parameter 

As a measure of the degree of intermolecular associa- 
tion, Hertz and co-workers [40, 41, 47] introduced a so- 
called "association parameter which is a weighted 
integral of the pair correlation function of the nuclei con- 
tributing to the dipolar relaxation process (in the present 
case nuclei in TBA-dl-solutions) and is defined as [48] 



A = 



174^2 



2 



(£)7°°(^%hhW4„^* 



(8) 



where a is the "closest approach distance of the interact- 
ing nuclei" and is usually assumed to be independent of 
the system's composition [48]. 

The correlation time T2 of the intermolecular dipole- 
dipole interaction is approximated by [48] 



^2, inter 



3D 



(9) 



where D is the self diffusion coefficient of the solute 
molecules. D can also be measured by NMR, using e.g. 
pulsed field gradient experiments [49]. 

Changes of the ^-parameter indicate short-range 
changes in the pair correlation function, which in the 
present study characterize the solvent-mediated interac- 
tion between the solute molecules. Enhanced associa- 
tion is identified by an increasing A-parameter: As the 
first neighbor peak of gunir) becomes sharper, the A- 
parameter increases 



TABLE I: Nonbonded interaction parameters used in the 
present study. The ^ refers to the ion parameter set of 
Heinzingcr [51], whereas the ^ refers to the parameters of 
Koncshan ct al. [52]. 
ing to aij = {au + crjj ; 



Lorentz-Berthelot mixing rules accord- 
,)/2 and 



: •/iiitjj were employed. 



Site 


?/|e| 


cr/nm 




OW 


-0.8476 


3.1656 


78.2 


HW 


-hO.4238 






CT 


-t-0.265 


3.50 


33.2 


CT(CH3) 


-0.180 


3.50 


33.2 


HC 


-fO.060 


2.50 


15.1 


OH 


-0.683 


3.12 


85.6 


HO 


-1-0.418 






Na^ 


+1.Q 


2.73 


43.06 


Cl^ 


-1.0 


4.86 


20.21 


Na^ 


-hl.O 


2.583 


50.32 




-1.0 


4.401 


50.32 



Using the definitions of Eq. 8 and Eq. 9, A is given in 
terms of NMR measurable quantities [48, 50] : 



A = 



1 



Ti 



, inter 



D 



(10) 



where pu is the number density of the ^H- nuclei in the 
system. Note that in our study pn is kept nearly con- 
stant, varying only the additional salt content, whereas 
in most NMR-studies [15-19, 48, 50] the concentration 
of the molecules, whose aggregation behavior is studied, 
is varied [18]. 

Besides detecting Ti^inter and D by NMR, we can also 
calculate these quantities, which are required for A, in- 
dependently from our simulations and determine the A- 
parameter for our model system exactly the same way 
as it is done from the experiment. Since the TBA-TBA 
pair-correlation fimction and the corresponding coordi- 
nation number is also available from MD, the simulations 
thus provide a "proof of concept" for a system behaving 
closely similar to the real system. 



C. MD simulation details 

We employ molecular dynamics (MD) simulations in 
the NPT ensemble using the Nose-Hoover thermostat 
[53, 54] and the Rahman- Parr inello barostat [55, 56] with 
couphng times tt = 1.5 ps and Tp = 2.5 ps (assuming the 
isothermal compressibility to be Xt = 4.5 10~^bar~^), 
respectively. The electrostatic interactions are treated 
in the "full potential" approach by the smooth particle 
mesh Ewald summation [57] with a real space cutoff of 
0.9 nm and a mesh spacing of approximately 0.12 nm and 
4th order interpolation. The Ewald convergence factor 
a was set to 3.38 nm~^ (corresponding to a relative ac- 
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TABLE II: Parameters characterizing the performed MD- 
simulation runs. All simulations were carried out at T = 298 K 
and P = 1 bar. The star * indicates the simulation run em- 
ploying the parameters of Koneslian et al. [52] for NaCl. 
For comparison a pure water simulation run of 500 SPCE 
molecules over 20 ns was performed yielding an average den- 
sity 998.4 kgm-^ 



7V(H20) 


1000 


1000 


1000 


1000 


N{TBA) 


20 


20 


20 


20 


iV(NaCl) 




10* 


10 


20 


Simul. length r/ns 


50 


100 


100 


100 


Density (p) / kg m~^ 


991.5 


1008.6 


1002.6 


1013.7 


(c(TBA))/moiri 


1.0170 


1.0044 


0.9984 


0.9810 


(c(NaCl)) /mol 




0.5022 


0.4992 


0.9810 



curacy of the Ewald sum of 10~^). A 2.0 fs timestep 
was used for all simulations and the solvent constraints 
were solved using the SETTLE procedure [58], while the 
SHAKE method was used to constrain the solute bond 
lengths [59]. All simulations reported here were carried 
out using the GROMACS 3.2 program [60, 61]. The 
MOSCITO suit of programs [62] was employed to gen- 
erate start configurations, topology files, and was used 
for the entire data analysis presented in this paper. Sta- 
tistical errors in the analysis were computed using the 
method of Flyvbjerg and Petersen [63]. For all reported 
systems initial equilibration runs of 1 ns length were per- 
formed using the Berendsen weak coupling scheme for 
pressure and temperature control (tt = rp = 0.5 ps) [64]. 



As in Refs. [25, 26] we study 0.02 mole fraction aque- 
ous solutions of t-butanol with and without presence of 
sodium chloride. The simulations were carried out for 
1 bar and 298 K. Our model system contains 1000 water 
molecules, represented by the three center SPCE model 
[65]. The flexible OPLS all-atom forcefield [66] is em- 
ployed for the 20 TB A- molecules. Here the bond lengths 
were kept fixed. Ten and twenty sodium chloride ion 
pairs were used to represent the salt solution. In order 
to check the influence of ion parameters variation, the 
two different parameter sets according to Hcinzinger [51] 
and Koneshan et al. [52] were employed. All non-bonded 
interaction parameters are summarized in Table I. To 
ensure proper sampling, and to allow an accurate deter- 
mination of the system's structural and dynamical prop- 
erties, the aqueous TBA solutions were studied for 50 ns, 
whereas the salt-solutions were monitored for 100 ns. The 
performed simulation runs and resulting concentrations 
are indicated in Table II. For comparison, a pure water 
system of 500 SPCE water molecules was simulated for 
20 ns for the same conditions. 



D. Experimental Details 

To experimentally determine the association behavior 
of TBA molecules in aqueous solutions we measured in- 
termolecular NMR relaxation rates, self diffusion coeffi- 
cients, and the densities of the aqueous solutions. All 
experimental data are summarized in Table HI. 

The observed relaxation rates depend on inter- as well 
as intra-molecular correlation functions. To extract the 
intermolecular rates, which are sensitive to the solute- 
solute association, we used the method of isotopic dilu- 
tion [67]. 

Since we are interested only in the hydrophobic methyl- 
protons, we deuterated the water and the hydroxyl group 
of the TBA. Isotopic dilution was performed by mix- 
ing (CD3)3COD (TBA-dlO) with (CH3)3COD (TBA- 
dl). We parameterize the dilution with the mole fraction 



[TBA-dl] 



[TBA-dl] + [TBA-dlO] ' 



The basic assumption of the isotopic dilution proce- 
dure is that the relaxation rate is given by the sum 
of an intramolecular term, which is independent of the 
dilution, and an intermoleciflar term, which is propor- 
tional to the concentration of the corresponding molecu- 
lar species. The contribution of the deuterated molecules 
can be taken as proportional to that of the protonated 
molecules, with a reduction factor [42] 

Ih{Ih + 1) 

The observed relaxation rate becomes therefore 

111 1 

+ ■ 



Ti Tin Ti 



+ 



-[{1 - a)xH + a]. (11) 



Here, Ti^intra denotes the intramolecular contribu- 
tion from protons within the methyl groups of the 
same molecule as the one being measured, Ti^inter 
the intermolecular contributions between diflferent TBA 
molecules, and Ti_o all other terms, such as paramagnetic 



TABLE 111: Experimental densities p, intermolecular relax- 
ation times Ti4ntor and self-diffusion coefficients D for TBA- 
dl in TBA-dl0/D2O/NaCl solutions. All experiments were 
carried out at T=298K at ambient pressure conditions. Also 
given are the TBA-dl concentrations and the obtained A- 
parameters. 



TBA-dl:D20:NaCl 2:100:0 2:100:1 



p/kgm 3 
c(TBA)/moir^ 

c(NaCl)/moir^ 

D/lO-V^s-i 
yl/10-3^m«s-2 



1083.9 
1.0059 

44.89 ± 1.74 
0.3962 
1.619 



1106.1 
0.9994 
0.5270 
41.02 ±4.73 
0.3818 
1.718 



2:100:2 
1122.5 
0.9882 
1.0403 
37.92 ± 0.13 
0.3792 
1.826 
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relaxation and interaction with other molecules such as 
D2O. To extract the intermolecular term, we measured 
the relaxation rate as a function of the isotopic dilution 
and fitted the measured data points to Eq. 11. 

The diffusion coefficients of TBA-dl were determined 
from the Pulsed Gradient Spin Echo (PGSE) experiments 
[49], where the gradient calibration was done using the 
diffusion coefficient of pure water [68] . 

The iH number density in the TBA-dl /D2O and TBA- 
dl/D20/NaCl solutions was obtained by measuring the 
mass density of the corresponding solutions with defined 
composition using a commercial Anton Paar oscillating 
U-tube density meter. 

The TBA-dl (99%) was purchased from Cambridge 
Isotope laboratories, the TBA-dlO (99%) from Isotec. 
The solvent D2O with the purity 99.96% was obtained 
from Merck KGaA. The solution was prepared by mea- 
suring the appropriate amount of each compound with a 
micropipette, and by weighing a corresponding amount 
of NaCl. The degassing process was done by the usual 
Freeze-Pump-Thaw technique, repeated several times un- 
til no gas bubbles didn't develop from the solution. Fi- 
nally, the samples were flame sealed. Relaxation and dif- 
fusion measurements were carried out at 600 MHz using 
a Varian Infinity Spectrometer system. All experiments 
were conducted under controlled temperature conditions 
at 25°C. 



III. RESULTS AND DISCUSSION 



A. Structural characterization of the aqueous TBA 
solutions 



The most prominent feature of the combined neutron 
scattering/EPSR work of Bowron and Finney is the ob- 
servation of a significant decrease of the height of the first 
peak of the central-carbon pair correlation function upon 
addition of sodium chloride [25, 26]. This decrease of the 
nearest neighbor peak is accompanied by an increase of 
the second neighbor peak located at a distance of about 
0.85 nm. From steric arguments and an analysis of their 
EPSR data, Bowron and Finney conclude that this pro- 
cess is according to the formation of chloride-bridged 
TBA-pairs. In Figure 1 we present the corresponding 
curves obtained from the present MD simulations. Sur- 
prisingly, exactly the opposite behavior is found. Upon 
addition of salt a notable increase of the first peak is ob- 
served, suggesting an enhanced TBA-TBA aggregation. 
Moreover, this increase is clearly more pronounced when 
increasing the salt concentration from about 0.5 molar to 
about 1 molar. Table IV provides a quantitative analy- 
sis of the pair correlation data in terms of coordination 
numbers 
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FIG. 1: Radial pair distribution functions between the TBA cen- 
tral carbon atoms (CC) in aqueous solutions at different salt con- 
centrations. The star * indicates the data belonging to parameter 
set of Koneshan et al [52] for NaCl. 




(12) 



FIG. 2: Several representative TBA-solvent pair atom-atom radial 
distribution functions obtained from the TBA:Water;NaCl 2:100:1 
solution. CC denotes the TBA central carbon atom, whereas CM 
and OH specify the methyl-carbon and hydroxyl oxygens, respec- 
tively. 



where pfj is the average number density of atom type 
(3. In order to provide comparability with the data ob- 
tained by Bowron and Finney, the values of r,jii,i and r^ax 
for the distance intervals given in the upper part of Ta- 
ble IV were taken from Ref. [26] . As the increasing peak 
height of the MD-data suggests, the TBA-TBA coordina- 
tion number increases with increasing salt concentration. 
Moreover, the increase is found to be significantly larger 
than the size of the associated errorbars. 

A more detailed picture of the structure of the aque- 
ous salt solutions of TBA is given in terms of selected 
solute/solvent pair correlation functions in Figure 2. To 
quantify the change of ion-solvation with increasing salt 
concentration, ion/water-oxygen and ion/TBA-oxygen 
(first shell) coordination numbers are given the in lower 
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TABLE IV; Coordination numbers for a 0.02 molar aqueous solution of t-butanol (TBA) with and without added sodium 
chloride. The coordination numbers are obtained by integrating over the distance interval indicated. For a direct comparison, 
the values of Ref. [26] were taken for rmin and rmax for the upper part of the table. For the data shown in the lower part of 
the table (first shell coordination numbers) integration over g{r) was performed to the first minimum of the ion-oxygen pair 
correlation function. CC-refers to the central carbon atom of TBA. OW denotes the water-yxygen while OH specifies the 
oxygen atom in the hydroxy 1- group. The star* indicates the parameter set for NaCl according to Koneshan et al. [52]. The 
neutron scattering data were taken from Ref. [26] . 



TBA:Water:NaCl 









2:100:0 




2:100:1 




2:100:2 


Atom-pair 


rmin/nm 


r-max/nm 


n-scatt. MD sim. 


n-scatt. 


MD sim. 


MD sim.* 


MD sim. 


CC-CC 


0.43 


0.75 


1.5 ±0.7 1.15 ±0.01 


0.8 ±0.5 


1.32 ±0.05 


1.27 ±0.01 


1.45 ±0.03 


CC-CC 


0.75 


1.00 


1.0 ±0.6 1.41 ±0.02 


1.8 ±0.8 


1.39 ±0.06 


1.39 ±0.02 
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OH-OH 


0.25 


0.35 
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Na-OW 


0.21 


0.30 




4.2 ± 1.1 


5.562 ±0.004 


5.582 ± 0.003 


5.454 ± 0.002 


Cl-OW 


0.28 


0.36 




4.9 ± 1.3 


5.715 ±0.003 


6.319 ±0.003 


5.629 ± 0.002 


OH-Na 


0.32 


0.55 




0.2 ±0.2 


0.098 ±0.001 


0.101 ±0.001 


0.191 ±0.002 


OH-Cl 


0.29 


0.38 




0.2 ±0.2 


0.0120 ±0.0002 


0.0134 ±0.0003 


0.0251 ±0.0004 


HO-Cl 


0.19 


0.32 
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FIG. 3: Schematic illustration of the two-dimensional cylin- 
drical pair distribution functions. Distribution of TBA and 
the ions around a central TBA molecule. The center of mass 
of the central TBA molecules is at {z — 0, r — 0) and the 
C-O-bond is aligned along the z-axis. 



part of Table IV. The changes in coordination numbers 
as the salt concentration doubles is only relatively small. 
However, apparently there is a certain tendency of the 
ions to slightly dehydrate, whereas the Na-OH coordina- 



tion number is almost unchanged and the Cl-OH coordi- 
nation slightly increases. 

In order to take the amphiphilic nature of TBA more 
properly into account wc also provide two-dimensional 
cylindrical pair correlation functions g{z,r), indicating 
the arrangement of molecules around a central TBA 
molecule using the notation of [46] 



1 



/Ntba Nfi 



NtbaPp 



J2 wi,) (13) 



i=i j=i 



where the z-axis is aligned along the TBA intramolecular 
unit-vector noH = ^oa/TOYi pointing from the center of 
mass to the hydroxyl oxygen. The index i runs over all 
TBA molecules, whereas the index j runs over a partic- 
ular subset of molecules (TBA, anions, or cations). The 
vector Tij = fj — fi represents the center of mass sepa- 
ration between particles i and j. The schematic shown 
in Figure 3 illustrates how the parameters r and z are 
defined. The 2D-distribution of TBA molecules around 
a central TBA molecule is shown in Figure 4. The ra- 
dial pair distribution functions in Figure 1 and the 2D- 
distributions in Figure 4 are interrelated and the radial 
distribution functions can be obtained by averaging over 
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FIG. 4: Two dimensional cylindrical center of mass pair distribution functions g(r, z) of TBA around a central TBA-molecule 
with the C-O- Vector pointing upwards, (a) 2 TBA : 100 Water, (b) 2 TBA : 100 Water : 1 NaCl; (c) Difference between the 
distribution functions shown in (b) and (a) (salt minus no salt). 
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FIG. 5: Two dimensional cylindrical pair distribution func- 
tions g{r, z) of the sodium (a) and chloride (b) ions around a 
central TBA-molecule with the C-O- Vector pointing upwards. 
From the simulatio with composition TBA:Water:NaCl of 
2:100:1. 



angles — arctan (r/z). Figure 4 reveals that the pre- 
peak in the CC-CC radial pair distribution function lo- 
cated a about 0.47 nm is due to hydrogen bonded TBA- 
TBA pairs. These pairs are identified by a separate dark 
spot at about r = 0.3 nm and z — 0.35 nm in close prox- 




FIG. 6: Schematic illustration of the two-dimensional cylin- 
drical pair distribution function g{r, z) of TBA molecules 
around a TBA — Cl~ contact pair with dTBA-ci- <0.64nm. 
The origin {z — 0,r =0) is located halfway between the cen- 
ters of mass of TBA and Cl~. The vector connecting TBA 
and Cl~ is aligned along the z-axis. 



imity to the hydroxyl-group in Figures 4a and 4b. As 
deduced from the radial distribution functions in Figure 
1, the addition of salt leads to an increased aggregation of 
TBA. The difference between the two-dimensional distri- 
bution functions with and without salt, shown in Figure 
4c, reveals that aggregation occurs predominantly on the 
methyl-group side of the TBA molecule: The region with 
negative z shows an increase in peak height of about 0.5, 
whereas on the hydrophilic side, the peak heights remain 
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0.2 0.4 0.6 0.8 1.0 1.2 
r/ nm 

FIG. 7: Cylindrical pair distribution function TBA around 
a contact pair formed by a TBA-molecule and a Cloride ion 
(TBA:bottom; Cl~: top). The geometry and the definition 
of r and z is indicated in Figure 6. 



almost unchanged. 

The distributions of the ions around a TBA-molecule 
are shown in Figure 5. Dark regions close to the hy- 
droxyl group indicate a significant stability of TBA-ion 
complexes. A remarkable difference is observed for ad- 
sorption of the different ions at the aliphatic side of TBA. 
The chloride ions tend to adsorb close to the Methyl- 
groups and are present as a grey shadow on the hy- 
drophobic side of the TBA molecule as shown in Fig- 
ure 5b, whereas the sodium completely tends to avoid 
that region. Consequently, there is practically no peak 
in the methyl-carbon-sodium radial pair correlation func- 
tion in Figure 2. This observation seems to be in line 
with the recently reported preferential anion-adsorption 
at the liquid-gas interface [69-71], which has structural 
similarities with the interface to a hydrophobic surface 
[11, 72, 73]. The tendency of the chloridc-ion to attach to 
the methyl groups might also explain the observed slight 
increase in the Cl-OH coordination number. A chloride- 
ion just gets more frequently attached to the OH-group 
since due to the preferential methyl-group-chloride inter- 
action the chloride sees a higher local TBA concentration. 

Finally, in order to prove whether there is a signif- 
icant amount of chloride-bridged TBA-TBA configura- 
tion we calculate the two-dimensional cylindrical (TBA — 
Cl~ ) . . . TBA correlation functions of a second TBA 
molecule around a TBA — CP contact pair with a center 



of mass distance less than 0.64 nm according to 
= -T-j 

iVTBA-Cl PTBA 

(JVtba-ci Wtba 
i=l 2 = 1 




Here fiTBA-ci is the unit vector describing the orienta- 
tion of a TBA — Cl~ contact pair. An illustration is given 
in Figure 6. As can be seen clearly from Figure 7, the vast 
amount of TBA molecules is located on the TBA-side of 
the TBA — CP -pair, practically ruling out a significant 
contribution of TBA - CP - TBA bridges. 

As the 2-dimensional cylindrical pair distribution func- 
tions indicate, the MD simulations furnish a scenario of 
an increased number of hydrophobic contacts in an aque- 
ous solution of TBA in the presence of sodium chloride, 
in contrast to the observations of Bowron and Finney. 
However, when taking a look at the published data, and 
considering the given size of errorbars, the n-scattering 
and MD simulation are not at all contradictory. A quan- 
titative comparison of the obtained coordination num- 
bers is given in Table IV. Almost all indicated values 
agree within the given errors, or are at least very close 
to each other. In fact, the large error in the n-scattering 
data for the central carbon (CC) coordination numbers 
for the first and second hydration shell, probably forbids 
a clear distinction between both scenarios just relying on 
the n-scattering data. 

Very recently Lee and van der Vegt have published a 
series of molecular dynamics simulations of TBA/ water 
mixtures with varying composition [74]. In their study, 
using the SPC model [75] for water, they found the orig- 
inal united-atom OPLS potential [76] for TBA less sat- 
isfactory and hence derived an improved set of model 
parameters for tertiary butanol, basically by adjusting 
the charges on the OH-group. Lee and van der Vegt re- 
port a too strong association of the TBA molecules for 
the original united-atom OPLS model. We would like to 
point out that reported effect might be also partly at- 
tributed to the use of the SPC water model, since an 
analogous solvent-model dependence has been recently 
observed for simple solutes (noble gases). Paschek [77] 
found for the SPC model a substantially stronger hy- 
drophobic association, which experiences also a signifi- 
cantly smaller temperature dependence. For our present 
study we choose the SPCE model on purpose, since it re- 
produces several water properties somewhat better than 
SPC, such as 00-radial distribution functions [78, 79], 
the self-diffusion coefficients [80], as well as the solvation 
entropy of small hydrophobic particles [77]. However, us- 
ing their models for an aqueous salt solution, preliminary 
data of Lee and van der Vegt indicate an association be- 
havior of TBA very much similar to our MD simulations: 
An enhanced hydrophobic TBA-TBA association upon 
addition of sodium chloride [81]. 
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TABLE V: Statistical analysis of the composition of the first solvation shell of the solvated Chloride-ion. n specifies the 
number of TBA-molecules in the solvation shell. P„ indicates the probability of finding a chloride ion with a solvation shell 
containing n TEA molecules. {Nw)„ indicates the average number of water molecules found in a solvation shell containing n 
TEA molecules. The left part of the tabic indicates data obtained directly from the MD-simulations. In the right part of the 
Table the Pj^"*' {Keq) indicates the combinatorial prediction of the probability of finding a chloride ion with a solvation shell 
containing n TEA molecules (see text for details). Here K^q is the equilibrium constant describing the a priori Water/TEA 
equilibrium for each hydration shell site. 
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B. A combinatorial picture of TB A/water-content 
equilibrium in the chloride solvation shell 

In this section we would like to elucidate the mech- 
anism that determines the composition of the chloride 
solvation shell by putting forward a simplified picture 
based on purely combinatorial arguments. This analysis 
will provide us also an estimate of the relative difference 
in binding free energy between TBA and water found 
in our MD simulation, as well as the one that would 
be necessary to create the Bowron-Finnney scenario of 
Chloride-bridged TBA-pairs. 

First of all we have statistically analyzed the struc- 
ture of the first solvation shell of a chloride ion as it is 
provided by our MD simulations. We have calculated 
the probability of finding a solvation shell with no TBA, 
one TBA molecule, two TBA molecules and so forth. 
As members of the first hydration shell are considered 
all TBA molecules with roH-ci < 0.39 nm and all wa- 
ter molecules with row-ci 0.39 nm (see Table IV and 
Figure 2). These data are shown in the left part of Table 
V. Given are results for MD simulations at two differ- 
ent salt concentrations. Moreover, we also provide the 
average number of water molecules (.^w)„ in a solvation 
shell accommodating n TBA molecules. 

The Cl-OW coordination number shown in Table IV, 
as well as the values for (iVw)^ given in Table V indicate 
that on average about seven to eight molecules (water 
and TBA) form the solvation shell of a chloride ion in 
the simulated solution. Now we would like to apply a 
simplified model of the solvation shell. Let us assume we 
have m sites (in our case: say seven) in the shell that 
can be either occupied by a water, or a TBA molecule. 
Then the probability that a site is occupied by a TBA 
molecule is according to 

-Keq ■'Vt + A/w 



where Nt and iVyy specify the number of TBA and water 
molecules in the simulation box, respectively. K^q repre- 
sents the constant describing the TBA/ water equilibrium 
per solvation shell site. In fact, Equation 15 represents 
the saturation limit for a Langmuir-type adsorption of a 
binary mixture (see Ref. [82] p. 80). Correspondingly, 
AG = —RT \nKcq represents the relative free energy 
change when replacing one water molecule by a TBA, 
and might be considered as the difference in "binding 
strength" . Given there are no further inter-site correla- 
tions within the solvation shell, the probability to find 
a chloride ion surrounded by exactly n TBA molecules 
follows a binomial distribution, and is given by 

P™('<1) (ifeq) = (1 - Pt)(™-") . (16) 

On the right hand side of Table V we illustrate three 
different cases: For the case of equal binding strength 
for water and TBA {K^q = 1), the probabilities P„ are 
just determined by the composition of the solution (the 
NT/N^J\f ratio). We would like to emphasize that even 
for this case the probability to find Cl-bridged dimers is 
below one percent. When comparing with data obtained 
directly from the MD simulation, we find them well repro- 
duced assuming Keq'vO.25. This indicates that a TBA- 
molccule is on average more weakly bound to a chloride- 
ion than a water molecule by about AG « 3.4kJmoP^. 
We would like to point out that due to the small value 
of pt ~ 5 X 10"^, the obtained probabilities P„ depend 
only very weakly on the variation of the number of sites 
m (say between 6 and 9) and are well approximated by 
a Poisson distribution. At last we would like to discuss 
a case, where the maximum in P„ is found for the case 
of Cl-bridged TBA-dimers (n = 2). To generate such 
a scenario, the equilibrium has to be shifted largely to 
the TBA-side with K^q ps 20, leading to a roughly 80 
times stronger binding of TBA compared to our MD- 
simulations, or an increase of about llkJmol^ on a 
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free energy seale. Coneluding, we would like to em- 
phasize that unless the ion-TBA binding is not substan- 
tially stronger than the ion-water interaetion, simply the 
vast majority of water molecules in the solution will pre- 
vent the formation of Cl-bridged TBA-dimers. Finally, 
although our employed model potential is certainly far 
from being perfect, there is actually no indication on 
what could cause such a strong specific chloride-HO inter- 
action. Moreover, the variation of the partial charges on 
the hydroxyl hydrogen (HO) varies only slightly among 
different forcefields: -|-0.418|e| for the all-atom OPLS 
model [66], +0.423|e| for van der Vegt's optimized united- 
atom OPLS model for TBA [74], and +0.435|e| for the 
united-atom TraPPE-forcefield by Chen et al. [83]. A 
specifically linearly chloride-bridged TBA pair would rep- 
resent a state of additional configurational order and 
hence even lower entropy. Favoring this particular con- 
figuration among all others would therefore require an 
increasingly lower energy for this type of complex. 
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FIG. 8: Log-log plot of the probability p{s) to find a TBA- 
cluster of size s, where s indicates the number of TBA- 

molcculcs in the cluster. The solid line indicates the behavior 
expected for random bond percolation on a SD-lattice close 
to the percolation transition [86]. 



C. TBA-clustering and "salting out" 

To determine whether the observed enhanced hy- 
drophobic aggregation upon addition of salt leads to a 
"salting out" of the TBA-molecules, we present a clus- 
ter size analysis of the TBA-aggregates. Such an ap- 
proach has been recently put forward to study the onset 
of a demixing transition for the case of binary water- 
tetrahydrofuran mixtures [84] and to structurally char- 
acterize aqueous methanol solutions [85]. Based on 
the TBA-TBA central carbon pair correlation function 
(shown in Figure 1) we consider two TBA molecules 
as "bound" when their central carbon C-C distance is 
smaller than 0.72 nm. We would like to point out that 
also hydrogen bonded TBA-TBA pairs are included by 
this definition. Based on that definition of intermolecular 
connectivity, we identify clusters of associated molecules. 
Figure 8 shows the log-log plot of the probability to find 
a cluster of a certain size s (number of TBA molecules) 
for the aqueous TBA solutions with increasing salt con- 
centration. First of all we would like to point out that 
for all shown distributions and for small cluster-sizes, the 
distribution approaches p(s) ~ s""^ with an exponent of 
r w 2.18, as observed for the case of random bond perco- 
lation on a 3D-lattice close to the percolation transition 
[86]. This seems to be well in accordance with the picture 
of an equilibrium of randomly associating and dissociat- 
ing clusters of molecules. For larger cluster sizes, how- 
ever, we observe a systematic deviation from the power- 
law behavior, which is due to the constraint of a finite 
number of TBA-molecules in our simulation. Please note 
that with increasing salt concentration an increased num- 
ber of large clusters appears. This is a direct evidence 
for a "salting out" effect upon increased salt concentra- 
tion. We would like to stress the observation, that also 
by visual inspection the solutions still appear to be ho- 
mogeneous and are apparently situated below a possible 
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FIG. 9: Distribution of the size of the largest clusters Smax 
(number of TBA molecules) per configuration as a function 
of salt concentration. 



phase separation. Moreover, for a system close to a 
phase separation, we would expect the appearance of a 
"hump" at large cluster-sizes. The relation of the build- 
up of such a narrow distribution of cluster-sizes showing a 
maximum close to the largest possible cluster size (in our 
simulation: 20 TBA molecules) with the onset of phase 
separation has been recently demonstrated for the case of 
tetrahydrofuran/water mixtures by Oleinikova et al.[84]. 
However, we cannot fully rule out that the tendency of 
the system to phase separate might be also suppressed by 
finite size effects. To further substantiate the apparent 
"salting out" tendency, we show in Figure 9 the distri- 
bution of cluster-sizes of the largest cluster detected in 
each configuration. The maximum of each of the distri- 
bution functions is located in the vicinity of 5 molecules. 
However, there is a clear shift towards larger cluster- 
sizes with increasing salt concentration. By promoting 
larger cluster-sizes with increasing NaCl concentration, 
the TBA-clustering analysis clearly indicates a "salting 
out" tendency of the aqueous TBA-solutions. 
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D. Self-difFusion coefficients 

We have determined the self diffusion coefficients for 

all particle types in the MD simulation from the mean 
square displacement according to 

where c represents the position of the molecules center of 
mass. In fact, the diffusion coefficients are obtained from 
the slope of the mean square displacement over the time- 
interval between 100 ps and 500 ps. The lower boundary 
has been chosen to be large compared to the the aver- 
age intermolecular association times. The diffusion coef- 
ficients obtained from the simulations arc given in Table 
VI. In addition, the experimental self-diffusion coeffi- 
cients for TBA-dl in in heavy water /salt solutions are 
shown in Table III. 

It is evident that the experimental diffusion coefficients 
of TB A are substantially smaller than the values obtained 
from MD simulations. This has to be largely attributed 
to the fact that heavy water is used as a solvent in the 
experiment, which has a significantly smaller diffusion 
coefficient {D = 1.768 x lO'^ s"^ at 298.25 K [87]) 
compared to H2O = 2.299 x IQ-'^ m^s^^ at 298.2 K 
[68, 88]). However, even when taking tliis effect into ac- 
count, the diffusion coefficient of TEA according to the 
MD simulations seems to be overestimated by about 40%, 
which could indicate an increased number of associated 
TEA- TEA pairs in the experiment. This discrepancy 
may at least partly also be attributed to a possible en- 
hanced "retardation effect" [89] of heavy water in the 
hydrophobic hydration shell of TEA. The "retardation 
effect" is clearly seen in the simulations as the slowing 
down of the water molecules in the TB A/water solution. 
However, given that the water retardation effect is related 
to the structuring of water in the hydrophobic hydration 
shell, it should rather be weaker in SPCE than in real 
water, since SPCE has been shown to underestimate the 
solvation entropy of hydrophobic particles [77] . 

In addition, we find that the decrease of the TEA 
diffusion coefficient upon addition of salt by 6.4% per 
moir^ NaCl is significantly smaller than by the 15.8% 
per moll~^ NaCl as obtained in the MD simulation. The 
overestimated salt-effect might be attributed to the po- 
tentially too small diffusion coefRcients observed for Na"*" 
and particularly Cl~. Scaling the ion self-diffusion co- 
efficients with the diffusion coefficient of water in the 
TE A/water system (as approximation), one would ex- 
pect diffusion coefficients of £'(Na''')R::il.l8 x 10~^ ms~^ 
and D{Cr) w 1.70 x lO'^ ms'^ for the c(NaCl) = 
1.0 moll~^ concentration (using diffusion coefRcients 
from Ref. [90] for the ions in aqueous solutions at the 
given concentration for 25°C). Although the interac- 
tion with the TEA molecules might also have a non- 
negligible effect, the MD simulations seem to underes- 
timate the self-diffusion coefficients of the ions by about 
13% and 46%, respectively. To confirm this observation 
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FIG. 10: ^-parameter rclatcd-quantitics as obtained from the 
MD-simulation: a) Inverse self-diffusion coefficient of TBA 
versus tfie intermolecular dipolar correlation time T2, inter 
of the aliphatic protons in TBA for the different MD simu- 
lations, b) TBA-TBA aliphatic carbon coordination number 
versus the j4-parameter. 

we have additionally calculated self-diffusion coefficients 
for a purely aqueous salt solution containing 500 SPCE 
molecules and 16 ion pairs at 298 K and 1 bar, with a con- 
centration of 1.67 mol. We obtain diffusion coefficients of 
£)(Na+) = 1.07X 10-9 ms-2, D{Cr) = 1.21 x 10"^ ms'^, 
and i:>(H20) = 2.18 x 10"*^ ms^^. 

The diffusion coefficients of the ions clearly indicate 
that there is a need for improvement of the employed ion 
parameter sets. Particularly the hydration-strength of 
the chloride ion seems to be overestimated by the present 
models. The overall agreement with the n-scattering data 
and the less satisfactory diffusion data, however, might 
just reflect the more pronounced sensitivity of kinetic 
quantities on details of the pair correlation function, such 
as the depth of the first minimum in the ion-water pair 
correlation function. This property is critically related 
to the water-exchange rate. 



E. TBA-TBA association and the A-parameter 

Given the uncertainties associated with the use of 
empirical forcefield-models, we have tried to find out, 
whether there is also experimental evidence for an en- 
hanced hydrophobic association scenario, as it is sug- 
gested by our MD simulations. Therefore we perform nu- 
clear magnetic relaxation experiments on aqueous solu- 
tions of tertiary butanol with varying salt concentration. 
Given a non-changing intermolecular dipolar correlation 
time r2, inter J an increasing relaxation rate of the methyl 
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TABLE VI: Self diffusion coefficients Aeif as obtained from the MD simulations for the TBA/Water- and TBA/Water/Salt- 
mbctures and for pure water. The star* indicates the parameter set for NaCl according to Koneshan et al. [52] . 

D/IQ-^ m^s-i 
TBA:Water:NaCl 



Particle 2:100:1* 2:100:2 2:100:1 2:100:0 0:100:0 



H2O 2.035 ±0.003 1.975 ±0.005 2.075 ± 0.005 2.197 ±0.006 2.642 ± 0.012 

TBA 0.751 ±0.009 0.707 ± 0.009 0.749 ±0.011 0.840 ± 0.015 

Na+ 0.987 ±0.018 1.041 ±0.013 1.091 ±0.017 

Cr 1.200 ±0.021 1.162 ±0.014 1.290 ± 0.020 



TABLE VII: Parameters characterizing the full H-H-dipolar correlation function (according to Eq.3) for the TBA aliphatic 
hydrogen nuclei as obtained from the MD-simulations. Given are both intra- and intermolccular contributions. The values 
indicate a by 20 % increased intermolccular rela^cation rate (Tj~^)inter in the presence of NaCl, whereas the intramolecular 
rate changes only by about 2 %. To obtain the correlation times r2 and prefactors (^'hh) = ^2(0) (see Eq.4), the correlation 
function was integrated numerically while the tail of the correlation function was fitted to a single exponential function (inter: 
between 120 and 200 ps; intra: between 30 and 50 ps). The corresponding intra- and intermolecular relaxation times were 
finally obtained from Eq.2. 
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T2, inter and the inverse self-diffusion coefficient of the so- 
lute molecules (Equation 9), which can be both obtained 
independently. To prove this we have calculated the full 
inter- and intramolecular dipole-dipole correlation func- 
tions according to Eq. 3 for the aliphatic hydrogens for 
all MD-simulations, assuming all other protons to be ex- 
changed by deuterons. Here only like-spin — ^ H dipo- 
lar interactions are considered. A quantitative descrip- 
tion of the intra- and intermolecular contributions as well 
as the calculated relaxation times Ti arc given in Table 
VII. The TBA self-diffusion coefficients for each of the 
simulated systems is given in Table VI. Figure 10a shows 
that the intermolccular correlation times r2, inter ^^nd the 
inverse self diffusion coefficient of TBA are indeed almost 
linearly related, establishing the A-paramctcr approach 
as a valid approximation in the present case. Moreover, 
as Figure 10b demonstrates, the ^-parameter from the 
MD-data is also almost linearly related to the TBA-TBA 
coordination number, provoking an interpretation of the 
A-parameter as an approximate direct measure of the 
TBA-TBA coordination number. 

The A-parametcrs obtained experimentally from the 
measured diffusion coefficients and relaxation rates are 
summarized in Table III. Both experimental and simu- 
lated A-Parametcrs arc shown as well in Figure 11. MD- 
simulation and experiment agree at least in a qualita- 
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FIG. 11: ^-parameter obtained from MD simulation and ex- 
periment as a function of salt concentration. 



hydrogen nuclei would directly indicate an association be- 
havior as observed in our MD-simulations. However, a 
decreasing relaxation rate would support the scenario ob- 
tained by Bowron and Finney. Since the correlation times 
arc likely to be changing, we follow the yl-paramcter ap- 
proach proposed by Hertz and co-workers, and discussed 
extensively in a previous section. The ^-parameter ap- 
proach is based on the assumption of a linear relation- 
ship between the intermolecular dipolar correlation time 
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tive sense: with increasing salt concentration the rising 
A-parameter indicates enhanced hydrophobic TBA-TBA 
aggregation. The quantitative difference of about 25% 
between the experimental and simulated A-parameters 
might be partially attributed to imprecisions of the em- 
ployed potential model, but might be based as well to a 
certain degree on the necessity for using D2O as a sol- 
vent in the NMR-measurements. Making use of the al- 
most linear relationship between the A-paramcter and 
the TBA-TBA coordination number, we would finally 
like to determine approximate TBA-TBA coordination 
numbers for the NMR experimental data. Doing so, 
we get coordination numbers of Nnn{TBA) — 0.975, 
1.040, and 1.111 for the solutions with salt concentra- 
tions O.0moir\ 0.5270moir\ and 1.0403 moir\ re- 
spectively. Although by about 25% smaller than the co- 
ordination numbers obtained from MD-simulation, and 
suggesting a smaller concentration variation, the experi- 
mentally obtained ^-parameters and coordination num- 
bers still tend to confirm an enhanced hydrophobic ag- 
gregation of the TBA molecules with increasing salt con- 
centration. Finally, we would like to point out that these 
"experimental" coordination numbers are also in agree- 
ment with the errorbars of the n-scattering data shown 
in Table IV. 



IV. CONCLUSIONS 

We have used a combination of molecular dynamics 
simulations and nuclear magnetic relaxation measure- 
ments to investigate the effect of salt (sodium chloride) 
on the association-behavior of tertiary butanol molecules 
in an aqueous solution. We have shown that the ap- 
plication of the so-called "A-parameter" -approach, pro- 
posed by Hertz and co-workers, employing solute-solute 
intermolecular ^ff -relaxation times and solute diffusion 
coefficients to determine the (relative) degree of associa- 
tion of solute molecules in aqueous solution, is well justi- 
fied in the present case. Moreover, our MD-simulations 
establish an almost linear relationship between the A- 
parameter and the TBA-TBA coordination number. 

Both MD-simulations and NMR-experiment tend to 
support a classical hydration and "salting out" picture of 



an enhanced tendency of forming hydrophobic contacts 
between the TBA- molecules in the presence of salt. An 
increasing salt concentration is hence found to strengthen 
the solute-solute hydrophobic interaction. Consequently, 
also the TBA-cluster-size distributions reveal a growing 
size of the TBA-aggregates as the salt concentration in- 
creases and therefore show directly the "salting out" ten- 
dency. Apparently, the TBA molecules behave closely 
similar to purely hydrophobic solutes, as recently shown 
by Ghosh et al.[8] for the case of hydrophobic methane 
particles dissolved in aqueous salt solutions. In the light 
of the scenario recently suggested by Koga and Widom 
[5, 10], the salt thus would lead to an increase of the 
excess chemical potential of the hydrophobic groups (re- 
duce their solubility) and would therefore provoke an en- 
hanced aggregation. Based on lattice model calculations 
Koga and Widom have proposed an almost inverse lin- 
ear relationship between the excess chemical potential 
of hydrophobic solutes and their hydrophobic interaction 
strength. 

Our detailed structural analysis of the simulation 
data does not provide any evidence for the presence of 
chloride-bridged butanol-pairs, as proposed by Bowron 
and Finney [25, 26]. Moreover, a combinatorial analy- 
sis of the composition of the chloride solvation shell re- 
veals that the formation of a significant amount of dimers 
is probably just prevented by the vast majority of wa- 
ter molecules in the solution. Finally, we would like 
to emphasize that although our results suggest a struc- 
turally completely distinct scenario, the molecular dy- 
namics simulations, as well as the coordination numbers 
obtained indirectly from our NMR-experimental data, 
appear to be largely within the experimental errorbars 
of the n-scattcring data and therefore seem to be com- 
patible with Bowron and Finney work [26]. 

Acknowledgments 

We acknowledge support from the Deutsche 
Forschungsgemeinschaft (FOR-436 and GK-298) and the 
University of Dortmund ( "Forschungsband Molekulare 
Aspekte der Biowissenschaften"). The authors are 
grateful to M. Holz for critically reading the manuscript. 



[1] C. Tanford, The Hydrophobic Effect: Formation of Mi- 
celles and Biological Membranes (John Wiley & Sons, 
New York, 1980), 2nd ed. 

[2] A. Ben-Nairn, Hydrophobic Interactions (Plenum Press, 
New York, 1980). 

[3] L. R. Pratt, Annu. Rev. Phys. Chem. 53, 409 (2002). 

[4] N. T. Southall, K. A. Dill, and A. D. J. Haymet, J. Phys. 
Chem. B 106, 521 (2002). 

[5] B. Widom, P. Bhimalapuram, and K. Koga, Phys. Chem. 
Chem. Phys. 5, 3085 (2003). 

[6] W. L. Masterton, D. Bolocofsky, and T. P. Lee, J. Phys. 



Chem. 75, 2809 (1971). 
[7] M. Kinoshita and F. Hirata, J. Chem. Phys. 106, 5202 

(1997). 

[8] T. Ghosh, A. E. Garcia, and S. Garde, J. Phys. Chem. 

B 107, 612 (2003). 
[9] T. Ghosh, A. Kalra, and S. Garde, J. Phys. Chem. B 
109, 642 (2005). 
[10] K. Koga, ,J. Chem. Phys. 121, 7304 (2004). 
[11] D. Huang and D. Chandler, Proc. Natl. Acad. Sci. USA 

97, 8324 (2000). 
[12] H. S. Ashbaugh and L. R. Pratt, physics/0307109v2. 



14 



[13] S. Rajamani, T. M. Truskctt, and S. Garde, Proc. Natl. 

Acad. Sci. USA 102, 9475 (2005). 
[14] D. Chandler, Nature (London) 437, 640 (2005). 
[15] M. Holz and M. Sorensen, Ber. Bunsenges. Phys. Chem. 

96, 1441 (1992). 
[16] M. Holz, R. Grunder, A. Sacco, and A. Meleleo, J. Chem. 

Soc. Faraday Trans. 89, 1215 (1993). 
[17] A. Sacco, F. M. De Cillis, and M. Holz, J. Chem. Soc. 

Faraday Trans. 94, 2089 (1998). 
[18] M. Mayele, M. Holz, and A. Sacco, Phys. Chem. Chem. 

Phys. 1, 4615 (1999). 
[19] M. Mayele and M. Holz, Phys. Chem. Chem. Phys. 2, 

2429 (2000). 

[20] G. W. Euliss and C. M. Sorensen, J. Chem. Phys. 80, 
4767 (1984). 

[21] D. T. Bowron, J. L. Finney, and A. K. Soper, J. Phys. 

Chem. B 102, 3551 (1998). 
[22] D. T. Bowron, A. K. Soper, and J. L. Finney, J. Chem. 

Phys. 114, 6203 (2001). 
[23] D. T. Bowron and J. L. Finney, Phys. Rev. Lett. 89, 

215508 (2002). 

[24] S. Dixit, J. Grain, W. C. K. Poon, J. L. Finney, and A. K. 

Soper, Nature (London) 416, 829 (2002). 
[25] D. T. Bowron and S. Di'az Moreno, J. Chem. Phys. 117, 

3753 (2002). 

[26] D. T. Bowron and J. L. Finney, J. Chem. Phys. 118, 
8357 (2003). 

[27] M. G. Cacace, E. M. Landau, and J. J. Ramsden, Q. Rev. 

Biophys. 30, 241 (1997). 
[28] F. Hofmeister, Arch. Exp. Path. Pharmakol. 24, 247 

(1888). 

[29] V. A. Parsogian, Nature (London) 378, 335 (1995). 
[30] B. Hribar, N. T. Southall, V. Vlachy, and K. A. Dill, J. 

Am. Chem. Soc. 124, 12302 (2002). 
[31] A. Geiger, Ber. Bunsenges. Phys. Chem. 85, 52 (1981). 
[32] R. Leberman and A. K. Soper, Nature (London) 378, 

364 (1995). 

[33] J. L. Finney and A. K. Soper, Chem. Soc. Rev. 23, 1 
(1994). 

[34] A. K. Soper, Chem. Phys. 202, 295 (1996). 

[35] A. K. Soper, Mol. Phys. 99, 1503 (2001). 

[36] D. T. Bowron, Philos. T. Roy. Soc. b 359, 1167 (2004). 

[37] J. L. Finney and D. T. Bowron, Curr. Opin. Colloid. In. 

9, 59 (2004). 

[38] J. L. Finney and D. T. Bowron, Philos. T. Roy. Soc. A 
363, 469 (2004). 

[39] M. G. Burke, R. Woscholski, and S. N. Yaliraki, Proc. 
Natl. Acad. Sci USA 100, 13928 (2003). 

[40] H. G. Hertz, Prog. NMR Spec. 3, 159 (1967). 

[41] H. G. Hertz and R. Tutsch, Ber. Bunsenges. Phys. Chem. 
80, 1268 (1976). 

[42] A. Abragam, The Principles of Nuclear Magnetism (Ox- 
ford University Press, 1961). 

[43] P. Wcstlund and R. Lynden-Bell, J. Magn. Reson. 72, 
522 (1987). 

[44] M. Odehus, A. Laaksonon, M. H. Levitt, and 
J. Kowalewski, J. Magn. Reson. Sor. A 105, 289 (1993). 

[45] E. O. Stejksal, D. E. Woessner, T. C. Farrar, and H. S. 
Gutowsky, J. Chem. Phys. 31, 55 (1959). 

[46] P. A. Egelstaff, An Introduction to the Liquid State (Ox- 
ford University Press, Oxford, 1992), 2nd ed. 

[47] A. L. Capparelli, H. G. Hertz, and R. Tutsch, J. Phys. 
Chem. 82, 2023 (1978). 



[48] K. J. Miiller and H. G. Hertz, J. Phys. Chem. 100, 1256 

(1996). 

[49] W. S. Price, Concepts Magn. Reson. 10, 197 (1998). 
[50] U. Wandle and H. G. Hertz, Z. Phys. Chem. 178, 217 
(1992). 

[51] K. Heinzinger, in Computer Modelling of Fluids Polymers 
and Solids, edited by C. R. A. Catlow, S. C. Parker, and 
M. P. Allen (Kluwer Academic Publishers, Dordrecht, 
1990), vol. C293 of NATO ASI Series, pp. 357-369. 

[52] S. Koneshan, J. C. Rasaiali, R. M. Lyndon-Bell, and S. H. 
Lee, J. Phys. Chem. B 102, 4193 (1998). 

[53] S. Nose, Mol. Phys. 52, 255 (1984). 

[54] W. G. Hoover, Phys. Rev. A 31, 1695 (1985). 

[55] M. Parrinello and A. Rahman, J. Appl. Phys. 52, 7182 
(1981). 

[56] S. Nose and M. L. Klein, Mol. Phys. 50, 1055 (1983). 
[57] U. Essmann, L. Perera, M. L. Berkowitz, T. A. Darden, 

H. Lee, and L. G. Pedersen, J. Chem. Phys. 103, 8577 

(1995). 

[58] S. Miyamoto and P. A. Kollman, J. Comp. Chem. 13, 

952 (1992). 

[59] J. P. Ryckaert, G. Ciccotti, and H. J. C. Berendsen, J. 

Comp. Phys. 23, 327 (1977). 
[60] E. Lindahl, B. Hess, and D. van der Spoel, J. Mol. Model. 

7, 306 (2001). 

[61] D. van der Spoel, E. Lindahl, B. Hess, A. R. van Buuren, 
E. Apol, P. J. Moulonhoff, D. P. Tioleman, A. L. T. M. 
Sijbers, K. A. Feenstra, R. van Drunon, et al., Gromacs 
User Manual version 3.2, www.gromacs.org (2004). 

[62] D. Paschek, MOSCITO4 molecular dynamics simulation 
package (2005), (www.moscitomd.de). 

[63] H. Flyvbjerg and H. G. Petersen, J. Chem. Phys. 91, 461 
(1989). 

[64] H. J. C. Berendsen, ,J. P. M. Postma, W. F. van Gun- 
steren, A. DiNola, and J. R. Haak, J. Chem. Phys. 81, 
3684 (1984). 

[65] H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, 

J. Phys. Chem. 91, 6269 (1987). 
[66] W. L. Jorgenson, D. S. Maxwell, and J. Tirado-Rives, J. 

Am. Chem. Soc. 118, 11225 (1996). 
[67] G. Bonera and A. Rigamonti, J. Chem. Phys. 42, 171 

(1965). 

[68] R. Mills, J. Chem. Phys. 77, 685 (1973). 

[69] E. M. Knipping, M. ,J. Lakin, K. L. Foster, P. Jung- 
wirth, D. J. Tobias, R. B. Gcrber, D. Dabdub, and B. J. 
Finlayson-Pitts, Science 288, 301 (2000). 

[70] P. Jungwirth and D. J. Tobias, J. Phys. Chem. B 104, 
7702 (2000). 

[71] P. Jungwirth and D. J. Tobias, J. Phys. Chem. B 105, 
10468 (2001). 

[72] H. A. Patcl, E. B. Naumann, and S. Garde, J. Chem. 

Phys. 119, 9199 (2003). 
[73] I. Brovchenko, A. Geiger, and A. Oleinikova, J. Phys. 

Condens. Matt. 16, S5345 (2004). 
[74] M. E. Lee and N. F. A. van der Vegt, J. Chem. Phys. 

122, 114509 (2005). 
[75] H. J. C. Berendsen, J. P. M. Postma, W. F. van Gun- 

steren, and J. Hermans, in Interrnolecular Forces, edited 

by B. PuUmann (D. Reidel, Dordrecht, 1981), pp. 331- 

338. 

[76] W. L. Jorgensen, J. Phys. Chem. 90, 1276 (1986). 
[77] D. Paschek, J. Chem. Phys. 120, 6674 (2004). 
[78] G. Hura, J. M. Sorenson, R. M. Glaeser, and T. Head- 
Gordon, J. Chem. Phys. 113, 9140 (2000). 



15 



[79] J. M. Sorenson, G. Hura, R. M. Glacscr, and T. Head- 
Gordon, J. Chem. Phys. 113, 9149 (2000). 

[80] D. van der Spoel, P. J. van Maarcn, and H. J. C. Berend- 
sen, J. Chem. Phys. 108, 10220 (1998). 

[81] N. F. A. van der Vegt (2005), private communication. 

[82] A. Ben-Naim, Statistical Thermodynamics for Chemists 
and Biochemists (Plenum Press, New York, 1992). 

[83] B. Chen, J. J. Potoff, and J. I. Siepmann, J. Phys. Chem. 
B 105, 3093 (2001). 

[84] A. Oleinikova, I. Brovchcnko, A. Geiger, and B. Guillot, 
J. Chem. Phys. 117, 3296 (2002). 

[85] L. Dougan, S. P. Bates, R. Hargreaves, J. P. Fox, 
J. Grain, J. L. Finney, V. Reat, and A. K. Soper, J. 



Chem. Phys. 121, 6465 (2004). 
[86] D. Stauffer and A. Aharony, Introduction to Percolation 

Theory (CRC Press, Boca Raton, Florida, 1994), 2nd ed. 
[87] W. S. Price, H. Ide, Y. Arata, and O. Soderman, J. Phys. 

Chem. B 104, 5874 (2000). 
[88] M. Holz, S. R. Hell, and A. Sacco, Phys. Chem. Chem. 

Phys. 2, 4740 (2000). 
[89] R. Haselmeier, M. Holz, W. Marbach, and H. We- 

ingartner, J. Phys. Chem. 99, 2243 (1995). 
[90] R. Mills and V. M. M. Lobo, Self- diffusion in electrolyte 

solutions, vol. 36 of physical sciences data (Elsevier, Am- 
sterdam, The Netherlands, 1989). 



